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ABSTRACT 



< 

Context. Two main scenarios for the formation of the Galactic bulge are invoked, the first one through gravitational collapse or 

<""j , hierarchical merging of subclumps, the second through secular evolution of the Galactic disc. 

Aims. We aim to constrain the formation of the Galactic bulge through studies of the correlation between kinematics and metallicities 
in Baade's Window (/ = 1°, b = -4°) and two other fields along the bulge minor axis (/ = 0°, b = -6° and b = -12°). 
Methods. We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral 
resolution of R=20000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body 

ffi ■ simulations of the Galactic bulge. 

Results. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. 2010: the metal-rich population 
presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc one. 
In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008), visible also in the kinematics, can 
be related to a varying mix of these two populations as one moves away from the Galactic plane, alleviating the apparent contradiction 
Q\ ' between the kinematic evidence of a bar and the existence of a metallicity gradient. 

Conclusions. We show evidences that the two main scenarios for the bulge formation co-exist within the Milky Way bulge. 
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S 1. Introduction lAthanassoul"3 (120051)1. After the bar forma t ion it heats in th e 



X 



3 d2QOj 
vertical direction (Combes & Sanders ( 198 1 ), Merrifield ( 1996)) 
Although the Milky Way bulge is our closest opportunity giving rise tQ the typical boxy/peanut aspect Observational 
to study in detail such a complex chemo-dynamical sys- data of individual stars in our Galaxy p rovide s evidence 
tem, its formation and evolution is still poorly understood. of both scenarios . The preS ence o f a bar in the inner 
Indeed the high extinction, the crowding, and the superpo- Galaxy has been first suggested by | de Vaucouleursl (l964) 
sition of multiple structures along the line of sight make from gas kinematics and confirme d since then by numer- 
studies of the inner Galactic regions challenging. Two main ous stu dies in cluding infrared luminosity distribution (e.g. 
scenarios have been invoked for bulge formation. The first B litz & Spergel (1991)), photometric data (e.g. Staneketal. 
one is gravitational collapse (Eggen et al. 1962) or hierar- (1994)> Babusiaux & Gilmore (2005), Benjamin et al. (2005)), 
chical merging of subclumps ^NogucluJ dl999j), | Aguerri et al. | OH/IR and Si0 -maser kinematics (e.g. Habing et al. (2006)), 



d2001J)). In this case the bulge formed before the disc and ste llar kinematics (e.g. fc angwala et al. (2009)) "and microlens 

the star-formation time-scale was very short. The resulting ing surveys ( _ Alcocket al. (2000)). Several triaxial struc- 

stars are old (> 10 Gyr) and have enhancements of a ele- tures may even coexis t withi n the Galactic inner reg ions 

ments relative to iron which are characteristics of classical (AJard (2 o6H), iNishivaml^n (I20051) . iLckTe z-Corredoii-a et al. 

bulges. The other scenario is secular evolution of the_disc ^of). The boxy aspect of the bulge, detected in near-infrared 

ftrovigh a bar forming a^seudo-bulge (see, e.g.|Combes et al. Hght profile (Dwek et al J [1995k, also ^ gues for a secular for . 

(11990|), i NormanetalJ dl996j), | Kormendy & Kenmcutt| (|2004f), mation of the bulge Qn the other hand, from medium and high 

resolution spectroscopic data enhancements of a elements have 

* Based on ESO- VLT observations 71.B-0617, 73.B-0074, and Paris been observed (M c William & Rich ( 1994), Zoccali et al. (2006), 
Observatory GTO 71.B-0196. 
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iFulbright etafl (l2007h . lLecureur et all (120071) ') suggesting a short 
formation time-scale. The discovery of a possible dwarf rem- 
nant in the bulge by Ferraro et al. (2009) also argues for a hier- 
archical merging scenario. A metallicity gradient has been ob- 
serve d in the minor-axis direction of the bulge (Zoccal i et al.l 
2008) favouring a bulge formation through dissipational col- 
lapse. Recently some observed similarities between the [a/Fe] 
bulge tr end and the t rend o f t he local thick stars {Melende z et alj 
(2008b, iRvde et all d2010l) . lAlves-Brito eTaTl OPJOb) suggest 
that the bulge and the local thick disc experienced similar 
chemical evolution histories. In what concerns the age, colour- 
magnitude diagrams show th at most of the Galactic bulge stars 
are old e r than lOGyr ( Ortolan i et al.ldl995l).lFeltzing & Gilmord 
( 120001) . IZoccali et al.l d2003l) . IClarkson et all d2008l)) a l thoug h 
an intermediat e-age population exists ( Ivan Loon et al. (2003), 
iGroeneweeen & Blommaert (2005)). 

Chemo-dynamical modelling of disc galaxy formation in 
a Cold Dark Matter (CDM) universe stresses the fact that 
both types of bulges can coexist in the same galaxy. In the 
Samland & Ge rhard! d2003l) simulation the bulge contains two 
stellar populations, an old population formed during the collapse 
pha se and a younger bar popul ation, differing in the [a/Fe] ra- 
tio. Nakasato & Nomoto (2003) suggest that the Galactic bulge 
may consist of two chemically different components: one rapidly 
formed through subgalactic clump merger in the proto-Galaxy, 
and the other one formed later in the inner region of the 
disc. The [Fe/H] abundance of the merger component tends 
to be smalle r than that of the second component. Recently 
Rahi mi et al.l {2010) simulated bulges formed through multiple 
mergers and analysed the chemical and dynamical properties of 
accreted stars with respect to locally formed stars: accreted stars 
tend to form in early epochs and have lower [Fe/H] and higher 
[Mg/Fe] ratios, as expected. 

Exploring correlations of abundances and kinematics in 
Baade's Window, Soto et al. (2007) suggest a transition in the 
kinematics of the bulge, from an isotropic oblate spheroid to 
a bar, at [Fe/H] = -0.5 dex. Their stu dy was based on 315 K 
and M giants with proper motions of fSpaenhau er et alj d 19921) . 
radial velocities of Sadler et al. ( 1996) and low-resolution abun- 
dances of Terndrurj et al. ( 1995) re-calibrated with the iron abun- 
dances scale of Fulbright et al.l (120 06). Similar conclusions were 
obtained earlier bv lZhao et alJ dl99 4) using a sample of 62 K gi- 
ants who pointed out the triaxiality of the bulge from kinematics 
in Baade's Window and the fact that the metal poor and metal 
rich populations were not drawn from the same distribution. 

In the present paper we analyse the correlations between 
kinematics and metallicity along three different minor-axis 
fields: Baade's Window ( b = -4°), b = -6° a nd b = -12°, for 
which [Fe /H] abundan ces dZoccali et al.l d2008l) hereafter Paper I 
and iHill et al.l (2010) hereafter Paper II) and radial velocities 
data for about 700 stars have been determined. Paper II obtained 
[Fe/H] and [Mg/H] metallicity distributions for red clump stars 
in Baade's Window and showed that the sample seems to be 
separated into two different populations in the metallicity dis- 
tributions. We here use the kinematic properties of the samples 
to confirm and constrain the nature of these two populations 
in Baade's Window and study their relative proportion change 
along the bulge minor axis fields of Paper I. The paper is or- 
ganized as follows: section 2 summarizes the data used in this 
study. In section 3 we analyse Baade's Window by combining 
our spectroscopic data with OGLE proper motion data. In sec- 
tion 4 we analyse the radial velocity versus metallicity trend 
along the bulge minor axis. In section 5 we compare our data 



with some published N-body models. Section 6 discusses our 
results. 



2. Data 

The different samples are described in detail in Paperl and 
in Paper II. Only a brief summary is given here. Paperl sam- 
ples consist of K giants (called RGB samples in what follows) 
observed in the following fields: Baade's Window (7 = l°.l, 
b = -4°.0: 194 stars); b = -6° (/ = 0°.2, b = -6°.0: 
201 stars) and b = -12° (/ = 0°.0, b = -12°.0: 99 stars). 
All the stars have photometric V, I data coming from differ- 
ent sources (for details see Paperl and Paper II) and 2MASS 
J,H,K data. Spectra have been obtained with VLT/FLAMES- 
GIRAFFE spectrograph in Medusa mode at resolution of about 
20 000. The sample of Paper II, observed in the same conditions, 
contains 219 red clump stars (called RC data in what follows) 
in Baade's Window (/ = 1.0°, b - -3°. 9). Iron abundances have 
been obtained with a mean accuracy of about 0.2 dex, depending 
on the value of [Fe/H], more metal-rich stars have the largest er- 
rors. In order to increase the statistics, we have combined the 
two samples of Baade's Window (RC sample of Paper II and 
RGB sample of Paperl). However due to a difference in the re- 
duction process, the two original MDs were not totally compat- 
ible. We therefore used in this combination the new reduction 
of the RGB sample made with the PapeiTI automatic reduction 
process (see section 5.1.4 of Paper II) which lead to two com- 
patible MDs which we can therefore combine in section [3] of 
this study . Howeve r we n ote that, according to the Besancon 
model (Robin et al. (2003), see figure 9 of Paper II and figure 
10 of Paperl), the two selections have slightly different mean 
distance and contamination. In section [4] the other minor axis 
fields at b — -6° and b - -12° are compared with the Baade's 
Window data and, to be consistent, we will use only the orig- 
inal RGB data of Baade's Window computed in Paperl in this 
comparison. 

Radial velocities have been obtained using the cross- 
correlation tool available in the GIRAFFE reduction pipeline 
(Rover et al. 2002). The individual spectra were cross-correlated 
with a box-shaped template corresponding to a K0 giant star. For 
each target, the barycentric radial velocities derived for each ex- 
posure are combined into an average velocity, and the standard 
deviation is used as an error estimate. For the different fields, 
the median error on the combined velocities range from 0.26 to 
0.43 km/s. Throughout this paper, V r stands for the heliocentric 
radial velocity. 

OGLE-II proper motions dSumi et alJ l2004) are only avail- 
able for the Baade's Window field. We removed stars with less 
than 100 data points used in the proper motion computation. 
Those removed stars are either near the edge of the CCD im- 
age or CCD defects or affected by blending. In the present paper 
the mas/yr has been converted in km/s assuming a distance to 
the Galactic Centre of 8 kpc (1 mas/yr = 37.9 km/s). The mean 
error on the proper motions of our sample is 1 .24 mas/yr, i.e. 47 
km/s. We worked on the relative values of the proper mo tions 
as provided in the catalogue. Section 6 of lSumi et al.l ([2004) de- 
tails how one can transform the catalogue into an inertial frame 
based on the measured proper motions for the Galactic Centre 
(GC). Our data cover fields BUL_SC45 and BUL_SC46, which 
do not s how any differenc e in the proper motions of the GC (ta- 
ble 3 of lSumi et al.l (2004)) so that the relative zero-points from 
these two fields do not need to be taken into account. 

Throughout this paper errors will be computed from a boot- 
strap analysis using 1000 samplings of the datasets. 
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3. Analysis of Baade's Window 

Our sample of 340 stars in Baade's Window is at present the 
largest one that has proper motions, radial velocities and homo- 
geneous iron abundances determinations. This sample combines 
both the Baade's Window dataset of Paper I (RGB sample) and 
the dataset of Paper II (RC sample) with metallicities derived in 
an homogeneous way by Paper II. We have co mputed the proper 
motion dispersions following the equations in Spaenhaue r et al.l 
(fl99l : 



°5 = r J -ni>-0 2 - 1 Z£ 

^ (« - 1) ^-f n *-J 



(1) 



with Hi and £,, being respectively the individual proper motions 
and their measurement errors. We assume that the measurement 
errors are uncorrected. For the radial velocity dispersion com- 
putation the errors can be neglected compared to the intrinsic 
dispersion. 

The covariance crf b , used in the computation of correlations 
between the different velocity components and in the the veloc- 
ity ellipsoid angles, also needs to be corrected from measure- 
ment errors due to the transformation from image coordinates to 
galactic proper motions. Indeed we have: 

f* I - laHa + lgf*6 

Hb = b ar i a +bsHs 

with the errors in the equatorial proper motions un-correlated 
(Sumi 2010, private communication), which implies: 



1 " 

<r% = covin,, n b ) ~ ~ ^J^bc,^ + kb s f si ) 

n /=i 



(2) 



with go, and £5, the individual measurement errors on the equa- 
torial proper motions. In Baade's Window l a b a - —l$bs - 
-0.426. 

Our Baade's Window radial velocity measurements 07 - 
111+4 km/s, (V r ) -9 + 6 km/s are in excellent agreement 
with previous studies: Howard et al. (2008): 07 — 1 12 + 10 km/s, 
(V r ) = -5 + 14 km/s, Rangwala et al. (2009): 07 = 112 + 3 
km/s, (V r ) = -1 + 5 km/s, and other measurements listed in 
Rangwala et alJ d2009h . The proper motion dispersions 07 = 
3.13 + 0.16 mas/yr and <Tb - 2.50 + 0.10 mas/yr are also in 
ex cellent agreemen t with the high accuracy HS T measurements 
of iKozIowski et alJ (120061) : 07 = 2.87 + 0.08, cr b = 2.59 ± 0.08 
mas/yr. 



3.1. kinematics versus metallicity 

Figure Q] plots the dispersion of velocity components along the 
line-of-sight (07), the galactic longitude (07) and the galactic 
latitude (crb) as a function of [Fe/H]. The figure shows that 07 
increases with metallicity, while 07 decreases. No significant 
variation of 07, is seen. To quantify if these variations are in- 
deed significant, we use the F-test to compare the variances of 
the velocity distributions for the first and last 33 % quantiles 
([Fe/H]< -0.14 and [Fe/H]> 0.30): 07. and 07 are indeed signifi- 
cantly different (p-value = 0.04 and 0.01 respectively), while the 
variation in crb is not si gnificant (p-valu e = 0.2). These results 
do not confirm those of lSoto et alJ (120071) . Their study, based on 
a compiled sample of proper motions, radial velocities and low- 
resolution abundances, showed only a rather shallow variation 
in <Tb, in the sense that the more metal poor stars have a higher 
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Fig. 1. Dispersion of the radial (07) and tangential velocities (07, 
07,) as a function of metallicity in Baade's Window by bins of 0.4 
dex. 



(Tb value. The discrepancies with the study of ISoto et al. (2007) 
may partly come from the sample itself. Our results are based 
on more accurate and homogeneous spectroscopic data. The dif- 
ferent contamination of the sample by foregrou nd objects may 
also b e another source of discrepancy. Recently, IClarkson et al.l 
(2008) from HST proper motions and estimated photometric 
parallaxes in the Sagittarius Window constructed the {l,b} ve- 
locity ellipse as a function of line-of-sight distance and demon- 
strated that its properties are sensitive to the distance to the con- 
sidered objects. They concluded that the depth-integrated veloc- 
ity ellipsoid of a small population of objects should be treated 
with caution. Ha lf of our sample are red clump stars while 
ISoto et aL (2007) studied M giants, so we expect our sample to 
be less biased toward closer stars. 

In what follows we analyse the shape and the orientation of 
the velocity ellipsoid for the data sets quoted in Table Q] From 
the whole sample, the observed anisotropy in o~ilo~b — 1.25 + 
0.09 is in good agr eement with th e results o f Spaenhauer et al.l 



(119921). iFeltzing & Johnsonl (120021) . iKuiiken & R ich (2002) and 
Kozlowsk i et al.1 (120061) . We detected also an anisotropy in 
o~b/o- r = 0.86 ± 0.05 and no significant one in <Ti/cr r = 1.07 ± 
0.07. Since the kinematic properties of the whole sample is an 
average of the properties of distinct populations, we analyse the 
shape of the velocity ellipsoid of the two metallicity samples. 

The metal poor sample ([Fe/H]< -0.14) shows a 07 higher 
than the other components, 07 =* 07,, and an anisotropy in 
ci/o~b - 1.34 + 0.17 (confirmed by a F-test with a 98% confi- 
dence) and in cr//cr r = 1 .39 + 0. 16 (a F-test cannot be performed 
here due to the large difference in the errors in the proper mo- 
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Table 1. Baade's Window velocity ellipsoid parameters. N is the number of stars, 07, 07, 07, are the velocity dispersions, C/ r , Cb r , 
Cib are the cross-correlation terms, /,, = 6>/ r , 8h r , Qm are the velocity ellipsoid angles. The values are computed for the full data set, 
the first and the last 33% quantiles of the metallicity distribution. The values predicted by the dynamical models of Zhao 1996 and 
Fux 1999, discussed in section 5, are also provided. The velocity ellipsoid parameters of the Fux 1999 model are provided for all 
the stellar particles as well as separately for the stellar particles originally in the disc and in the old spheroid. 
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Fig. 2. V r and/i/ distributions in Baade's Window. 



tion a nd radial velocity measurements). Following IZhao et all 
(1996), the observed anisotropics for the metal poor sample 
(07/07, — (TilcTr) may be interpreted in terms of rotation broad- 
ening: as Baade's Window passes close to the minor axis of the 
bulge, only 07 should be broadened by the rotation, which is in- 
deed what is observed. However, a possible contamination with 
foreground thick disc stars higher than expected would also con- 
tribute to an increase in 07. 

For the metal rich sample ([Fe/H]> 0.30) the shape of the 
velocity ellipsoid is significantly different from the metal poor 
sample. We observe for the metal rich population 07 > 07 > 07,, 
and that 07. is higher (with 96% confidence) and 07 is lower (with 
99% confidence) than the observed values for the metal poor 
sample. In the case of a bar, the velocity dispersion along the 
bar major axis is expected to be much larger than its azimuthal 
dispersion which is itself lower than its vertical dispersion. If the 
bar was pointing nearly end-on, that would mean 07. > 07 > 07, 
which is what we observe. The anisotropy in <Tblcr r — 0.78 + 
0.08 observed for the metal rich sample may be explained by the 
flattening of the bar driven population. 
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Fig. 3. Vertex deviation /,, in Baade's Window as a function of 
metallicity by bins of 0.4 dex. 



The Baade's Window velocity ellipsoid shows the character- 
istics of a non-axisymmetric system. We now analyse the orien- 
tation of the velocity ellipsoid. We first examine the correlations 
between the different velocity components. The obtained cross- 
correlation terms (C xy = cr^l{(r x cr y ), cr 2 xy being the covariance) 
are given in Table [TJ 

Cib is higher than the results of Kozlowski e t al.l (I2006I) 
(Cu, = -0.10 + 0.04) and Rattenb urv et all d2007l) (Cm = 
-0.16 + 0.03). We note that Rattenburv et al. (2007) used the 
same OGLE-II data but selected only stars with errors in the 
proper motion smaller than 1 mas/yr, introducing a bias towards 
closer stars but reducing the impact of the correction described 
in equation [2] 

A decrease with metallicity of the correlations between the 
velocity components seems to be present in our sample. We now 
quantify the observed correlations in terms of velocity ellipsoid 
angles in Table[TJwith: 



tan(2<9„) 



2ai 



(3) 



In the following we will call the vertex deviation l v - 0; r , 
which is the orientation of the axis of the velocity dispersion 
tensor in the radial-longitudinal velocity plane illustrated in Fig. 
[2]for the whole sample. Table[TJand Figure[3]show that the metal 
poor sample is consistent with no vertex deviation (Pearson's 
correlation p-value = 0.1) while the metal rich sample shows 
a significant vertex deviation (/,, = -36° + 5 with a Pearson's 
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correlation p- val ue of 5 x 10~ 5 ). This confirms the results of 
ISoto et all(l2007l) and is a clear indication that the metal-rich part 
of our sample is kinematically under the influence of the bar. The 
metal-rich sample may consist essentially of stars formed by the 
secular disc evolution. Bars are intrinsically non-isotropic, the 
intrinsic velocity ellipsoid is longer in the radial direction (e.g. 
Zhao et all (1 19961) 1 in agreement with our results. Figure [3] indi- 
cates that the bar influence can already be seen at [Fe/H] 0.1 

dex. But we note that a small fraction of stars showing a correla- 
tion have a large influence o n the overal l correlation of a sample. 
We used the particles of the iFuxl (11999) model described in sec- 
tion 5 to test that indeed less than 30% of stars of disc/bar parti- 
cles are needed to introduce a significant correlation (p-value < 
0.05) in the disc/bar plus spheroid sample. 

Our derived v alue of 6/b = —22 ° + 6° is consistent with 
the results of ' darkson' et alj d2008h (6fe = -34° + 8°) and 
iKozlowski et al . (2006) (6ib — -24°). It shows a significant vari- 
ation with metallicity. The Pearson's correlation test indicates a 
correlation with 99% confidence for the metal rich sample while 
the metal poor sample correlation is not significant (p-value = 
0.1). 

We measure for the first time 8i, r . The Pearson's correlation 
test does not show a significant correlation between fib and V r 
for the metal poor sample (p-value=0.3), while it finds the cor- 
relation to be significant for the metal rich sample with 99.9% 
confidence. 

Following Paper II, we applied the SEM MUL Gaussian 
components decomposition algorithm (iCeleux & DiebollH l986) 
to the [Fe/H] distribution of the full Baade's Window sample 
and find a decomposition fully compatible with the Paper II 
results (which were based on the Red Clump sample only): 
the decomposition identifies clearly two populations, in roughly 
equal proportions, a metal-poor component (46+3% of the sam- 
ple, < [Fe/H] > =-0.32+0.02) with a large dispersion (0.23+0.01 
dex) and a metal-rich component (54±3% of the sample, 
<[Fe/H]>=0.30±0.01) with a small dispersion (0.11+0.01 dex). 
The variation of the kinematic properties between the metal rich 
and the metal poor parts of the sample studied here suggests dif- 
ferent formation scenarii for those populations. The metal poor 
component of our sample shows 07 higher than o> and 07,, o~r — 
o- b , a strong anisotropy in 07/0-j and 07 /o> and no correlations 
between the velocity components. This population can be inter- 
preted as a spheroidal population with the velocity anisotropy 
due mostly to rotation that broadens the observed 07 component. 
The metal rich population shows significant anisotropy and cor- 
relation between all the velocity components. In particular the 
strong vertex deviation indicates that it can be interpreted as a 
bar-driven population. 

3.2. Kinematics versus distance 

We used our red clump stars sample to probe the kinemat- 
ics against distance. Red clump stars should indeed allow 
to measure the velocity shift between the two bar streams 
dMao & Paczvhskil (l2002l) . lRangwala et al.1 (120091) 1. As the stars 
form, the bar streams in the same sense as the Galactic rotation 
and the bar being in the first quadrant, the stars on the near side 
of the bar are expected to go towards us, while stars in the far 
side should move away from us. The actual velocity shifts be- 
tween these two streams constrains the bar orientation angle. 

The Besancon model confirms that our red clump sample 
should be nicely centred around the Galactic Centre (see figure 
9 of Paper II). The magnitude / is less sensitive than magnitude 
K to both the metallicity and the age for ages older than ~4 Gyr 
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Fig. 4. Distribution of the / magnitude and V r in the Baade's 
Window Red Clump sample. The dotted lines represent the first 
and last 25% quantiles of the / magnitude distribution. The filled 
dots corresponds to stars with [Fe/H]> -0.09 used in Table|2] 



dSalaris & Girardil20 02). Moreover the / magnitude is more sen- 
sitive to the extinction which itself increases with distance, lead- 
ing the / magnitude to be much more sensitive to distance than 
K. No d ifferential e xtinction pattern is seen in the region of our 
sample (ISumill2004l) . 

We have cut our sample into bright and faint stars on the first 
and last 25% quantile of our / magnitude distribution (Fig. |4| 
and computed the mean velocities. Without selection in metal- 
licity, the Welch Two Sample t-test does not indicate that the 
mean values between the bright and the faint samples are signif- 
icant enough (p-value = 0.1 for V r , 0.1 for fi/ and 0.8 for fib). By 
removing the metal poor stars (selecting [Fe/H]> -0.09, which 
corresponds to the 66% quantile of the [Fe/H] distribution) the 
Welch Two Sample t-test becomes significant on the radial ve- 
locities (p-value = 0.02 for V r , 0.2 for fii and 0.4 for fib)- The 
resulting values are presented in Table [2] The numbers go in the 
expected direction which is a negative V r and a positive fi\ for 
the near stream and a positive V r with a negative fii for the far 
stream. Note that the mean proper motions quoted in Table|2]are 
relative. By removing the metal rich stars from the sample (se- 
lecting [Fe/H]<0.33, which corresponds to the 66% quantile of 
the [Fe/H] distribution) the Welch Two Sample t-test does not 
distinguish differences between the faint and the bright sample 
(p-value = 0.5 for V r , 0.2 for fi/ and 0.9 for fib). This again con- 
firms our interpretation of the metal-rich population of Baade's 
Window being associated to the bar and the metal-poor popula- 
tion to a spheroidal population. 

This shift was n ot seen in the radial velocities of 
Rangwala et al.l (12009!) . They argue that this may be due to 
Baade's Window population being the superposition of an old 
spheroidal population and a bar population, which we here con- 
firm. 
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Table 2. Baade's Window Red Clump mean velocities for the faint and bright stars in magnitude / (first and the last 25% quantiles) 
for stars with [Fe/H] > -0.09 (33% quantile of the [Fe/H] distribution). 



N 



(km/s) 



(mas/yr) 



(mas/yr) 



all 134 

/ < 15.05 37 
/ > 15.26 30 



7: 

-32: 

38: 



10 

20 

22 



0.4: 

0.7: 

-0.2: 



0.3 
0.5 
0.6 



-0.1. 

-0.4: 
0.2: 



0.2 
0.4 

0.5 



For the RGB sample, no correlation of the kinematics with 
the magnitude / can be seen. We note that the RGB sample is 
biased towards closest stars. According to the Besan§on model, 
the median distance difference between the RC and the RGB 
sample is expected to be 1 kpc (see figure 9 of Paper II for the 
RC and figure 10 of Paperl for the RGB sample). This bias in 
distance is not significant in the mean kinematics of the RGB 
sample (V r - 3 + 8 km/s, m - 0.6 + 0.3 mas/yr, /j. b = 0.5 + 0.2 
mas/yr). 

4. Radial velocity versus metallicity along the 
bulge minor axis 

We are now wondering whether the two distinct stellar popula- 
tions found in the previous section are present all along the bulge 
minor axes. We analyse the radial velocity distribution behaviour 
of the fields of Paper I, presented in Table [3] To be consistent 
with the b - -6° and b - -12° metallicity distribution, we will 
use only the RGB sample of Baade's Window with the origi- 
nal metallicity distribution function of Paperl. Paperl showed a 
variation of the metallicity distribution with height z above the 
Galactic mid-plane (Fig. |5}. In this section we combine metal- 
licity data with radial velocity data. Figure [6] shows the radial 
velocity dispersion as a function of metallicity. 

We first compare the kinematic behaviour along the bulge 
minor axis (TableO with previous results from the literature. We 
saw in section 3 that our Baade's Window radial velocity disper- 
sion is in excellent agreement with previous measurements. At 
b - -6° our radial velocity disper sion is lo wer than the BRAVA 
measure (o> = 108+7 km/s, Howard et al. (2008)). The decrease 
of the radial velocity dispersion with galactic latit ude we observe 
is con sistent with the SiO maser measurements of Izumiura et al. 
(119951) . 

All the fields show a negligible skew, while the kurto- 
sis are d ifferent. This variati on of the kurtosis was not de- 
tected in iHoward et all d2008l) . At b - -12° the distribution 
is significantly pointy (with 99.3% confidence according to the 
Anscombe-Glynn kurtosis test), indicating that the kinematics 
are significantly affected by the disc. At b — -6° the kurtosis be- 
comes consistent with zero. In Baade's Window the distribution 
is significantly flattened (with 93% confidence in the RGB sam- 
ple and 99.8% in the red -clump sample). Sharpies et al. ( 1990) 
and lRangwala et all d2009l) also measured a skew consistent with 
zero and kurtosis significantly negative in Baade's Window, in- 
dicating that the distributions are flat-topped rather than peaked. 
Rangwala et al. (2009) concludes that this seems to be consistent 
with a model of the bar with stars in elongated orbits forming 
two streams at different mean radial velocities, broadening and 
flattening the total distribution. 

We analyse now kinematic data versus metallicity. The mean 
radial velocities do not show significant variation with metal- 
licity in any of our fields. Figure [6] shows that the velocity 
dispersion at the rich end decreases significantly with latitude. 
However the velocity dispersion on the metal-poor end does not 



vary significantly with latitude. Summing all stars with [Fe/H] 
< -0.5 dex in all fields leads to cr r — 97 ± 7 km/s. 

The field b - -6° is close enough to Plaut's Window (/ = 0°, 
b - -8°) so that we can compa re our results to the proper mo- 
tion study of IVieira et al.l d2007l) . The Besangon model indicates 
a mean distance of our selec t ed bul ge stars to be 6.5 kpc, com- 
patible with the IVieira et al.l (2007) selection. Their metallicity 
distribution is compatible with ours. They do not find neither any 
change in the proper-motion dispersion as a function of metal- 
licity. They found 07 = 102 + 3 km/s and 07, = 88 + 3 km/s 
at b — -8°, we obtain o> = 83 km/s at b — -6°, leading to a 



full picture coherent with the predictions of the model of Zhao 
(1996) of 7, ~ o> and 7 > 07, at these latitudes (see the bottom 
ofFig.6of|Zhag|<[l996|)). 

We have seen in the previous section that Baade's Window 
kinematic behaviour as a function of metallicity can be inter- 
preted as a mix of two populations, a metal poor component with 
kinematics that can be associated to an old spheroid population 
and a metal rich component with bar-driven kinematics. Under 
this light we can interpret the variation with galactic latitude of 
both the metallicity distribution function (Fig. [5j and the kine- 
matic as a function of metallicity (Fig. [6) as the bar population 
disappearing while moving away from the plane. At high lat- 
itudes the foreground disc component dominates the metal rich 
part of the kinematic behaviour. The metal poor component asso- 
ciated with the old spheroid stays present along the bulge minor 
axis. 

We have obtained an estimation of the variation of the pop- 
ulations with latitude by using the SEMMUL Gaussian com- 
ponents decomposition in the metallicity distribution of the 
Paperl samples. Table|4]gives the decomposition of the Baade's 
Window RGB sample. The proportion of the metal-rich com- 
ponent is higher than in the red clump sample of PapeiTI and its 
spread in metallicity is higher also, most probably due to a differ- 
ence in the sample selection which is biased towards stars closer 
to the Sun in the RGB sample. Although both samples are not 
exactly on the same scale at the high metallicity end (cf section 
5. 1 .4 of Paper II), the mean metallicities for the metal-poor pop- 
ulation is extremely similar in both samples (see also the decom- 
position of the full Baade's Window sample on the same metal- 
licity scale at the end of section 3.1). At b — -6° the Wilks' test 
allows to keep a solution with 3 components presented in Table 
|5]rather than the 2 components one. Population A and popula- 
tion B could correspond to the population A and B observed in 
Baade's Window. The mean metallicities are coherent although 
their spread is smaller. The radial velocity dispersions of popula- 
tion A are identical. The radial velocity dispersion of population 
B decreases at b — -6° as expected for a bar-like kinematic 
behaviour (see Fig. [8] and associated text). Population C repre- 
sents only 6 + 2% of the sample with a low mean metallicity of 
—1.1 + 0.1 dex and a high velocity dispersion of 127+26 km/s, 
which could therefore be associated to the halo. The Besancon 
model prediction of only 0.4 ± 0. 1 % of halo star in the sample 
could therefore have been underestimated. This population could 
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Table 3. Summary of the bulge minor axis fields (Paper I) characteristics, zee 1S the height of the field along the bulge minor axis, 
assuming a distance Sun - Galactic Centre of 8 kpc. N is the number of stars. (V r ) is the mean heliocentric velocity. 



Field 


/ 
C) 


b 

C) 


ZGC 

(pc) 


N 


(Vr) 

(km/s) 


(km/s) 


skew 


kurtosis 


b = -4" 
b = -6° 
6 = -12° 


1.1 

0.2 
0.0 


-4.0 

-6.0 

-12.0 


606 

844 

1700 


194 

201 

99 


5±7 
-10±6 
-14 ±8 


104 ±5 
83 ±4 
80 ±8 


-0.07 ±0.1 
0.01 ±0.1 
0.38 ± 0.4 


-0.54 ± 0.2 

-0.19 ±0.2 

1.71 ±0.7 



U 



-12° 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

[Fe/H] 



- b = -4° 


I 


I 


i 

I J 


i 


- b = -6° 








- 


I 


I 


I 


- b = -12° 










II 


I 










J 



-1.5 -1.0 -0.5 0.0 0.5 

[Fe/H] 



1.0 



Fig. 5. Distribution of the metallicity for the different galactic lat- 
itudes of Paper I 



Fig. 6. Dispersion of the radial velocity for the different galactic 
latitudes as a function of metallicity by bins of 0.4 dex. 



also have been hidden in the population A at b — -4°. Selecting 
all stars with [Fe/H]< -0.9 in our 3 fields we obtain 16 stars 
with a radial velocity dispersion of o> = 116 + 21km/s, which 
is coherent with the solar neighbourho od velocity dispersion s 
measured in this metallicity range (e.g. IChiba & Beer s (2000): 
a~u ~ 110+10 km/s) containing both thick disc and halo stars. 
We are not in a position to clearly associate this population ei- 
ther to the halo or to a metal-poor thick disc in our inner galactic 
samples. SEMMUL did not converge on the b - -12° field due 
to the smaller number of stars and the higher contamination with 
thin, thick discs and halo stars expected in this field. At b - - 12° 
the metal rich component present alb — -4° and b - -6° seems 
to have fully disappeared. The metal rich velocity part of the ve- 
locity dispersion corresponds to a disc like component (see next 
section) while the metal poor part shows a velocity dispersion 
still coherent with the metal poor population of b — -4° and 
b - -6°, although we cannot distinguish a spheroid and a thick 
disc contribution. 



5. Comparison with dynamical models 

In this section we combine the full Baade's Window sample of 
section 3 with the b=-6° and b=-12° fields of section 4. 

We first compared our kinematics with the model of Zhao 
(Il996l) in Baade's Win dow (T able [B and along the minor axis 
(Fig. 0). The model of lZhaol (fl996) is a 3D steady-state stel- 
lar model using a generalized Schwarzschild technique consist- 
ing of orbital building blocks within a bar and a disc poten- 
tial. We find a very good agreement between our data and this 
model. Figure [7] shows also the very good agreement of our ra- 
dial velocity distribution along the minor axis with the BRAVA 
data dHoward et alj d2008l) . IHoward et ail d2009t) ). These com- 
parisons are done on global kinematics while we have shown 
here that several populations are present in these data. 

The N-body dynamical model of Fux ( 1999) allows to com- 
pare the kinematic properties of the particles that were originally 
in the disc to those being in the old spheroid. This model is a 
3D self-consistent, symmetry-free N-body and smooth particle 
hydrodynamics code. It contains 3.8 million particles: a dark 
halo, a disc and a spheroid and a gas component. Following 
IHoward et alj d2009l) we have used the Fux model cl0t2066 
which assumes a bar angle of 20° and a distance to the Sun 
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Table 4. SEMMUL Gaussian components decomposition of field b - -4°, RGB sample only. The radial velocity dispersion <x r is 
given in km/s. 



Pop [Fe/H] 



C[Fe/H] 



% 



&r 



A -0.31 ±0.05 0.39 ±0.02 33 ± 3 
B 0.13 ±0.02 0.23 ±0.01 67 ± 3 



109 ±7 



Table 5. SEMMUL Gaussian components decomposition of field b — -6°. The radial velocity dispersion cr r is given in km/s. 



Pop [Fe/H] 



&[FelH\ 



% 



cr r 



C -1.09 ±0.07 0.24 ±0.01 6±2 127 ± 26 

A -0.27 ±0.02 0.24 ±0.01 64 ± 3 83 ± 5 

B 0.14 ±0.02 0.13 ±0.01 30 ± 3 70 ± 6 



This work 
Howard et al. 2008 
Howard ct al. 2009 
Zhao 1996 




Fig. 7. Radial velocity dispersion along the bulg e mino r axis 
compared to the BRAVA data (iHoward et all (120081) and 
iHoward et al.1 ((2009)) and the model of ZhaoJ ((199J) 
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Fig. 8. Radial velocity dispe rsion along the bulge minor axis 
compared with the model of iFuxl d!999l) . The metal poor and 
metal rich population are defined by the 33% and 66% quantile 
of the [Fe/H] distribution. 



of 8 kpc. We have selected particles at a distance from the Sun 
6 < r < 10 kpc in a cone selection of 0.5 degree radius for 
\b\ < 7° and in a cone selection of 1 degree for \b\ > 7° to allow 
to work with a sample large enough away from the plane (this 
larger cone leads to 209 particles at b - -12°). We computed 
statistical uncertainties of the model by bootstrap of 2-3 km/s in 
the velocity dispersions and 5° in the velocity ellipsoid angles 
for the disc/bar, 20° for the spheroid. 



Table[T]shows the predictions of the model of Fux ( 1999) in 
Baade's Window. In this field the model contains 50% of stars 
from the spheroid and 50% from the disc/bar population, equiv- 
alent to the SEMMUL Gaussian components decomposition we 
obtained for the red-clump sample in Paper II. The velocity dis- 
persions of the model are slightly higher than our observed ones 
and the velocity components correlations are lower. The veloc- 
ity dispersions of the metal rich population corresponds to the 
disc/bar particles except for <x/ which is higher in the model. 
Both the metal rich observations and the disc/bar particles show 
a correlation of the velocity components but the correlations 
in the observations are higher. The vertex deviation is however 
identical in the model and in the observations. The metal poor 
component show lower velocity dispersions in cr r and crt, than 
the spheroid model particles. As expected the model does not 



show any correlations in the velocity components of the spheroid 
particles, in line with our metal poor sample. 

Figure [8] shows the strong decrease in radial velocity dis- 
persion predicted by the Fux model for the disc/bar population 
while moving away from the plane. In the plane the predicted 
dispersion is larger for the disc/bar than the spheroid one, due 
to the influence of the bar driven orbits. At high latitude the cr r 
dispersion is small for the bar/disc component, following a disc 
like behaviour. The spheroid population keeps the same velocity 
dispersion along the minor axis. In Fig. [8] we have over-plotted 
the Fux model radial velocity dispersions with our data for the 
full sample, the metal rich and the metal poor parts, as defined 
by the 33% and 66% quantiles in the [Fe/H] distribution. For 
the metal rich part we observe a strong decrease of the metallic- 
ity dispersion while going away from the plane, coherent with 
the behaviour of the disc/bar particles of Fux model. For the 
metal poor part we observe a constant velocity dispersion, as 
for the spheroid particles of the Fux model, but with a mean ve- 
locity dispersion of about 20 km/s lower than the model. It is 
this difference of velocity dispersion for the metal poor compo- 
nent which leads to a too high global velocity dispersion of the 
Fux m odel compared to the BRAVA data, leading Howard et al. 
(2009) to an interpretation different from ours: at b - -8° the 
radial velocity distribution is compatible with the disc/bar com- 
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ponent of the Fux model without the need for an old spheroid 
component. However with a velocity dispersion for the spheroid 
of ~ 100 km/s instead of ~ 120 km/s the global velocity disper- 
sion would be coherent with the full Fux model (spheroid+disc). 
The use of the Besancon model adds support to the interpretation 
th at the decre ase in cr r with latitude, seen in disc/bar particles 
of lFuxl {1999) and in our metal-rich samples, can be due to the 
disc replacing the bar in the sample: according to the Besancon 
model, the thin disc contamination of our sample at b — -6° was 
expected to be of 10% with a cr r - 62 km/s and 19% at b = -12° 
with cr r = 47 km/s. 

6. Discussion 

Our analysis of the kinematics as a function of metallicity in 
Baade's Window shows that our sample can be decomposed in 
two distinct populations for which we suggest different forma- 
tion scenarii. The metal poor component does not show any cor- 
relation between its velocity components and is therefore consis- 
tent with an isotropic rotating population. PapeiTI showed that 
this component is enriched in [Mg/Fe]. We interpret this popu- 
lation as an old spheroid with a rapid time-scale formation. The 
metal rich component shows a vertex deviation consistent with 
that expected from tags a population with orbits supporting a 
bar. Paper II showed that this component has a [Mg/Fe] near so- 
lar. We interpret this population as a pseudo-bulge formed over 
a long time scale through disc secular evolution under the ac- 
tion of a bar. This pseudo-bulge is gradually disappearing when 
moving away from the Galactic plane. 

In this context we can give a new consistent interpretation of 
the metallicity gradient in the bulge. A metallicity gradient is in- 
deed visible in the bulge whe n observing further away from the 
plane than Baade's Window dFrogel et al] d!999h and Paper I), 
while in the inner regi ons (\b\ < 4°) n o gr adient in m etallicity 
has been found dRamfrez et alJ d2000h and iRich et alJ d2007l) ). 
This can be understood if the bar as well as the old spheroid 
population are both present in the inner regions, leading to a 
constant metallicity at \b\ < 4°, while the bar influence gradu- 
ally fades further away from the plane than Baade's Window. At 
high latitudes only the old spheroid r emains: the mean metal- 
licity of the outer bulge measured by Ibat a & Gilmord {1995) 
of [Fe/H]~ -0.3 dex corresponds very well to our metal poor 
population. This scenario is also consistent with the distribu- 
tion of bulge globular clusters of Valenti et al. (2009) who found 
no evidence for a metallicity gradient but all their clusters with 
[Fe/H]> -0.5 dex are located within (\b\ < 5°). 

In what concerns the age of the two populations, we ex- 
pect the spheroid component to be old while the pseudo-bulge 
component may contain both the old stars of the inner disc 
redistributed by the bar and younger stars whose formation 
has been triggered by the bar gas flow. Ivan Loon et al.l d2003l) 
found that although the bulk of the bulge population is old, 
a fraction of the stars are of int ermediate age (1 to 7 Gyr). 
Groenewegen & Blommaert (2005) observed Mira stars of ages 
1-3 Gyr at all lat i tudes from -1.2 to -5.8 in the OGLE-II data. 
lUttenthaler et alJ d2007l) found 4 bulge stars with ages lower than 
3 Gyr at a latitude of b = -10°. IBensbv etalj d2009l) found 
3 microlensed bulge dwarfs with ages lower than 5 Gyr. 87% 
of the variable stars detected by Kouzuma & Yamaokal {2009) 
are distributed within \b\ < 5° and most of them should be 
large-amplitude and long-period variables such as Mira vari- 
ables or OH/IR stars. This intermediate age population has 
been shown to t race t he Galactic bar dvan Loon et al. ( 2003), 
llzumiura et al] d 19951) . iGroenewegen & Blommaert (2005), 



Kouzuma & Yamaoka (2009)), although providing a larger bar 
angle (~ 40°) than studies based on older tracers such as red 
clump stars (~ 20°, Stanek et al. ( 1994), Babusiaux & Gilmore 
d2005h ). This discrepancy in the bar angle could well be ex- 
plained if the old tracers probed a mix of spheroid and bar struc- 
tures while the young tracers only probe the bar one (although 
biases on the longitude area surveyed needs also to be taken into 
account, see Nishivama et al. (2005)). If this intermediate age 
population was associated to a part of the bar component, their 
presence in the CMDs would decrease while going away from 
the plane as the main bar co mponent and w ould therefore be a 
small fraction of the CMP of IZoccali et~aH d2003l) at b = -6°. 
IClarkson et al.1 (2008) obtained a proper motion decontaminated 
CMD with a well defined old turn-off in an inner field (I — 1°, 
b - -3°). However we would expect an intermediate age pop- 
ulation associated with the bar to be metal rich, which, due to 
the age-metallicity degeneracy, would imply that thi s po pula- 
tion could be hidden in the CMD of Cla rkson et al] d2008l) if 
its contribution is small enough compared to the bulk of the 
bulge population. The new filter combination proposed by the 
ACS Bulge Treasury Pr ogramme to break the age-metallicity- 
temperature degeneracy (Bro wn et al.l l2009) should provide new 
lights on this issue. 

We note that in Baade's Window, neither the kinematics nor 
the chemistry allows to distinguish what we call the old spheroid 
to the thick disc. The mean metallicity of the solar neighbour- 
hood thick disc is however lower (e.g. |Fuhrmann (2008J) de- 
rived [Fe/H] = -0.6 and [Mg/H] = -0.2) than the mean metal- 
licity of our metal poor p opulation (fFe/H l = -0 - 27 dex and 
[Mg/H1 = -0.04, Paper II). [ Melendez etal] d200 8). Rvde et al.l 
(2010J), IBensbv et all d2009l) and lAlves-Brito et al. (2010} ob- 
served similarities between the metallicity of the bulge and the 
meta llicity of thic k disc stars for metal poor stars. The sample of 
Rvde et all d2010h contains 1 1 of our stars all with less than so- 
lar metallicity. Simulations of the formation of thick stellar discs 
by rapid interna l evolution in unstable, gas-rich, clumpy discs 
(Bournaud et al. 2009) show that thick discs and classical bulges 
form together in a time-scale shorter than 1 Gyr which allow to 
explain the observed abundance similarities. 

The coexistence of classi cal and pseudo bulg e has been ob- 
served in external galaxies dPrugniel et al.l d200ll) . iPeletier et al.l 
(2007), Erwin (2008)) and obtained by N-bod y simulations 
(Samland & Gerhard (2003), Athanassoula (2005)). The chem- 
ical and dynamical model of Nakasato & Nomoto (2003) sug- 
gests the presence in the bulge of two chemically different com- 
ponents as we found, one formed quickly through the subgalac- 
tic clump merger in the proto-Galaxy, and the other has formed 
gradually in the inner disc. But they do not have a bar in their 
model. They fitted well the kinematics of Minniti (1996) who 
observed at / = 8°, b - 7° a decrease of <x r with metallicity, 
which i s coherent w ith what we obs erved at b = -12° and the 
fact that [Nakasato & Nomoto (200J) define d the bulge radius as 
R<2kpc. The chemo-dynamical model of ISamland & Gerhard! 
(2003) predicts the different characteristics of our sample: their 
total bulge population contains two stellar populations: a metal 
rich population ([Fe/H]>0.17) with [a/Fe]<0 associated with the 
bar, and an old population that formed during the proto-galactic 
collapse, with a high [a/Fe] and a [Fe/H] corresponding to the 
"thick disc" component. Their model also predicts the resulting 
apparent metallicity gradient along the bulge minor axis (their 
Fig. 13). 

Our study highlights the importance to combine metallicity 
to 3D-kinematic information to disentangle the different bulge 
populations. This approach needs to be extended to various 
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galactic longitudes. Gaia will not only provide those but also 
allow to determine the distances (probing the different structures 
along the line of sight and removing distance induced biases) 
and work on an impressively large sample of un-contaminated 
bulge stars. 

Acknowledgements. MZ and DM are supported by FONDAP Center for 
Astrophysics 15010003, the BASAL Center for Astrophysics and Associated 
Technologies PFB-06, the FONDECYT 1085278 and 1090213, and the 
MIDEPLAN MilkyWay Millennium Nucleus P07-021-F. 



References 

Aguerri, J. A. L., Balcells, M., & Peletier, R. F. 2001, A&A, 367, 428 
Alard, C. 2001, A&A, 379, L44 

Alcock, C, Allsman, R. A., Alves, D. R., et al. 2000, ApJ, 541, 734 
Alves-Brito, A., Melendez, J., Asplund, M., Ramirez, I., & Yong, D. 2010, ArXiv 

e-prints, 1001.2521 
Athanassoula, E. 2005, MNRAS, 358, 1477 
Babusiaux, C. & Gilmore, G. 2005, MNRAS, 358, 1309 
Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2005, ApJ, 630, L149 
Bensby, T., Feltzing, S., Johnson, J. A., et al. 2009, ArXiv e-prints, 0911.5076 
Blitz, L. & Spergel, D. N. 1991, ApJ, 379, 631 
Bournaud, E, Elmegreen, B. G., & Martig, M. 2009, ApJ, 707, LI 
Brown, T. M., Sahu, K., Zoccali, M., et al. 2009, AJ, 137, 3172 
Celeux, G & Diebolt, J. 1986, Rev. Stat. Appl., 34, 35 
Chiba, M. & Beers, T. C. 2000, AJ, 1 19, 2843 
Clarkson, W., Sahu, K., Anderson, J., et al. 2008, ApJ, 684, 1110 
Combes, E, Debbasch, E, Friedli, D., & Pfenniger, D. 1990, A&A, 233, 82 
Combes, E & Sanders, R. H. 1981, A&A, 96, 164 
de Vaucouleurs, G 1964, in IAU Symp. 20: The Galaxy and the Magellanic 

Clouds, 195-+ 
Dwek, E., Arendt, R. G, Hauser, M. G, et al. 1995, ApJ, 445, 716 
Eggen, O. J., Lynden-Bell, D., & Sandage, A. R. 1962, ApJ, 136, 748 
Erwin, P. 2008, in IAU Symposium, Vol. 245, IAU Symposium, ed. M. Bureau, 

E. Athanassoula, & B. Barbuy, 113-116 
Feltzing, S. & Gilmore, G. 2000, A&A, 355, 949 
Feltzing, S. & Johnson, R. A. 2002, A&A, 385, 67 

Ferraro, E R., Dalessandro, E., Mucciarelli, A., et al. 2009, Nature, 462, 483 
Frogel, J. A., Tiede, G P., & Kuchinski, L. E. 1999, AJ, 117, 2296 
Fuhrmann, K. 2008, MNRAS, 384, 173 

Fulbright, J. P., McWilliam, A., & Rich, R. M. 2006, ApJ, 636, 821 
Fulbright, J. P., McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1 152 
Fux, R. 1999, A&A, 345, 787 

Groenewegen, M. A. T. & Blommaert, J. A. D. L. 2005, A&A, 443, 143 
Habing, H. J., Sevenster, M. N, Messineo, M., van de Ven, G, & Kuijken, K. 

2006, A&A, 458, 151 
Hill, V., Lecureur, A., Gomez, A., et al. 2010, submitted, (Paperll) 
Howard, C. D., Rich, R. M., Clarkson, W., et al. 2009, ApJ, 702, L153 
Howard, C. D., Rich, R. M., Reitzel, D. B., et al. 2008, ApJ, 688, 1060 
Ibata, R. A. & Gilmore, G E 1995, MNRAS, 275, 605 
Izumiura, H, Deguchi, S., Hashimoto, O., et al. 1995, ApJ, 453, 837 
Kormendy, J. & Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603 
Kouzuma, S. & Yamaoka, H. 2009, AJ, 138, 1508 
Kozlowski, S., Wozniak, P. R., Mao, S., et al. 2006, MNRAS, 370, 435 
Kuijken, K. & Rich, R. M. 2002, AJ, 124, 2054 
Lecureur, A., Hill, V., Zoccali, M., et al. 2007, A&A, 465, 799 
Lopez-Corredoira, M., Cabrera-Lavers, A., Mahoney, T. J., et al. 2007, AJ, 133, 

154 
Mao, S. & Paczyriski, B. 2002, MNRAS, 337, 895 
McWilliam, A. & Rich, R. M. 1994, ApJS, 91, 749 
Melendez, J., Asplund, M., Alves-Brito, A., et al. 2008, A&A, 484, L21 
Merrifield, R. M. 1996, in Barred Galaxies, IAU Colloquium N157, ed. R. Buta, 

D. A. Crocker, B. C. Elmegreen, ASP, San Francisco, p. 179 
Minniti, D. 1996, ApJ, 459, 579 
Nakasato, N. & Nomoto, K. 2003, ApJ, 588, 842 
Nishiyama, S., Nagata, T, Baba, D., et al. 2005, ApJ, 621, L105 
Noguchi, M. 1999, ApJ, 514, 77 

Norman, C. A., Sellwood, J. A., & Hasan, H. 1996, ApJ, 462, 114 
Ortolani, S., Renzini, A., Gilmozzi, R., et al. 1995, Nature, 377, 701 
Peletier, R. E, Falcon-Barroso, J., Bacon, R., et al. 2007, MNRAS, 379, 445 
Prugniel, P., Maubon, G, & Simien, E 2001, A&A, 366, 68 
Rahimi, A., Kawata, D., Brook, C. B., & Gibson, B. K. 2010, MNRAS, 401, 

1826 
Ramirez, S. V., Stephens, A. W., Frogel, J. A., & DePoy, D. L. 2000, AJ, 120, 

833 



Rangwala, N., Williams, T. B., & Stanek, K. Z. 2009, ApJ, 691, 1387 
Rattenbury, N. J., Mao, S., Debattista, V. P., et al. 2007, MNRAS, 378, 1165 
Rich, R. M., Origlia, L., & Valenti, E. 2007, ApJ, 665, LI 19 
Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2003, A&A, 409, 523 
Royer, E, Blecha, A., North, P., et al. 2002, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 4847, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. J.- 
L. Starck & E D. Murtagh, 184-194 
Ryde, N., Gustafsson, B., Edvardsson, B., et al. 2010, A&A, 509, A260000+ 
Sadler, E. M., Rich, R. M., & Terndrup, D. M. 1996, AJ, 1 12, 171 
Salaris, M. & Girardi, L. 2002, MNRAS, 337, 332 
Samland, M. & Gerhard, O. E. 2003, A&A, 399, 961 
Sharpies, R., Walker, A., & Cropper, M. 1990, MNRAS, 246, 54 
Soto, M., Rich, R. M., & Kuijken, K. 2007, ApJ, 665, L31 
Spaenhauer, A., Jones, B. E, & Whitford, A. E. 1992, AJ, 103, 297 
Stanek, K. Z., Mateo, M., Udalski, A., et al. 1994, ApJ, 429, L73 
Sumi, T. 2004, MNRAS, 349, 193 

Sumi, T, Wu, X., Udalski, A., et al. 2004, MNRAS, 348, 1439 
Terndrup, D. M., Sadler, E. M., & Rich, R. M. 1995, AJ, 1 10, 1774 
Uttenthaler, S., Hron, J., Lebzelter, T, et al. 2007, A&A, 463, 251 
Valenti, E., Ferraro, E R., & Origlia, L. 2009, MNRAS, 1935 
van Loon, J. T., Gilmore, G E, Omont, A., et al. 2003, MNRAS, 338, 857 
Vieira, K., Casetti-Dinescu, D. I., Mendez, R. A., et al. 2007, AJ, 134, 1432 
Zhao, H, Rich, R. M., & Biello, J. 1996, ApJ, 470, 506 
Zhao, H, Spergel, D. N., & Rich, R. M. 1994, AJ, 108, 2154 
Zhao, H. S. 1996, MNRAS, 283, 149 

Zoccali, M., Hill, V, Lecureur, A., et al. 2008, A&A, 486, 177, (Paperl) 
Zoccali, M., Lecureur, A., Barbuy, B., et al. 2006, A&A, 457, LI 
Zoccali, M., Renzini, A., Ortolani, S., et al. 2003, A&A, 399, 931 



